PART 3 

INTRODUCTION TO ENGINEERING HEAT TRANSFER 



Introduction to Engineering Heat Transfer 



These notes provide an introduction to engineering heat transfer. Heat transfer processes set hmits 
to the performance of aerospace components and systems and the subject is one of an enormous 
range of application. The notes are intended to describe the three types of heat transfer and provide 
basic tools to enable the readers to estimate the magnitude of heat transfer rates in realistic aerospace 
applications. There are also a number of excellent texts on the subject; some accessible references 
which expand the discussion in the notes are listen in the bibliography. 
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1.0 Heat Transfer Modes 



Heat transfer processes are classified into three types. The first is conduction, which is defined 
as transfer of heat occurring through intervening matter without bulk motion of the matter. Figure 
1.1 shows the process pictorially. A solid (a block of metal, say) has one surface at a high 
temperature and one at a lower temperature. This type of heat conduction can occur, for example, 
through a turbine blade in a jet engine. The outside surface, which is exposed to gases from the 
combustor, is at a higher temperature than the inside surface, which has cooling air next to it. The 
level of the wall temperature is critical for a turbine blade. 



Thigh 




Heat "flows" to right ( q) 




Solid 



Figure 1.1: Conduction heat transfer 



The second heat transfer process is convection, or heat transfer due to a flowing fluid. The 
fluid can be a gas or a liquid; both have applications in aerospace technology. In convection heat 
transfer, the heat is moved through bulk transfer of a non-uniform temperature fluid. 

The third process is radiation or transmission of energy through space without the necessary 
presence of matter. Radiation is the only method for heat transfer in space. Radiation can be 
important even in situations in which there is an intervening medium; a familiar example is the heat 
transfer from a glowing piece of metal or from a fire. 

tAucfcfy fypinf^ 

How do we quantify the contribution of each mode of heat transfer in a given situation? 
(MP HT.l) 

2.0 Conduction Heat Transfer 

We will start by examining conduction heat transfer. We must first determine how to relate the 
heat transfer to other properties (either mechanical, thermal, or geometrical). The answer to this is 
rooted in experiment, but it can be motivated by considering heat flow along a "bar" between two 

heat reservoirs at Ta, Tb as shown in Figure 2.1. It is plausible that the heat transfer rate g , is a 
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function of the temperature of the two reservoirs, the bar geometry and the bar properties. (Are there 
other factors that should be considered? If so, what?). This can be expressed as 

Q = fi (Ta , Tb , bar geometry, bar properties) (2. 1) 

It also seems reasonable to postulate that Q should depend on the temperature difference Ta - Tb. If 
Ta - Tb is zero, then the heat transfer should also be zero. The temperature dependence can therefore 
be expressed as 

Q [ (T'a - Tb), Ta, bar geometry, bar properties] (2.2) 




Figure 2.1: Heat transfer along a bar 



An argument for the general form of f2 can be made from physical considerations. One 
requirement, as said, is f2= if Ta= Tb. Using a MacLaurin series expansion, as follows: 



f(^T) = f(0) + 



df 



d(AT) ^ 

If we define AT= Ta - Tb and f=f2, we find that (for small Ta - Tb), 



(2.3) 



f2(TA-TB) = Q = f2(0) + 



df2 



^(Ta-Tb) 



[Ta-Tb)- 



(2.4) 



Ta-Tb=0 



We know that f 2(6) = . The derivative evaluated at Ta = Tb (thermal equilibrium) is a measurable 

df. 

property of the bar. In addition, we know that Q>0 ifT^> Tg or —7 — — — r > 0. It also seems 

^(Ta - Tb) 

reasonable that if we had two bars of the same area, we would have twice the heat transfer, so that 
we can postulate that Q is proportional to the area. Finally, although the argument is by no means 
rigorous, experience leads us to believe that as L increases Q should get smaller. All of these lead 
to the generalization (made by Fourier in 1807) that, for the bar, the derivative in equation (2.4) has 
the form 
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3/2 



= f. (2.5) 

Ta-Tb=0 



In equation (2.5), A: is a proportionality factor that is a function of the material and the 
temperature, A is the cross-sectional area and L is the length of the bar. In the limit for any 
temperature difference AT across a length Ax as both L, Ta-Tb ^ 0, we can say 

e = yt^&l^=-fc4&^ = -M^. (2.6) 
L L dx 

A more useful quantity to work with is the heat transfer per unit area, defined as 

4 = ^ • (2-7) 

A 

The quantity q is called the heat flux and its units are Watts/m^. The expression in (2.6) can 
be written in terms of heat flux as 

q = -k^. (2.8) 
dx 

Equation 2.8 is the one-dimensional form of Fourier's law of heat conduction. The 
proportionality constant k is called the thermal conductivity. Its units are W / m-K. Thermal 
conductivity is a well-tabulated property for a large number of materials. Some values for familiar 
materials are given in Table 1 ; others can be found in the references. The thermal conductivity is a 
function of temperature and the values shown in Table 1 are for room temperature. 



Table 2.1: Thermal conductivity at room temperature for some metals and non-metals 



Metals 


Ag 


Cu 


Al 


Fe 


Steel 


k [W/m-K] 


420 


390 


200 


70 


50 



Non-metals 


H2O 


Air 


Engine oil 


H2 


Brick 


Wood 


Cork 


k [W/m-K] 


0.6 


0.026 


0.15 


0.18 


0.4 -0 .5 


0.2 


0.04 
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2.1 Steady-State One-Dimensional Conduction 



Insulated 
(no heat transfer) 

I I 

dx 



-> X 



Figure 2.2: One-dimensional heat conduction 



For one-dimensional heat conduction (temperature depending on one variable only), we can 
devise a basic description of the process. The first law in control volume form (steady flow energy 
equation) with no shaft work and no mass flow reduces to the statement that 2g for all surfaces = 
(no heat transfer on top or bottom of flgure 2.2). From equation (2.8), the heat transfer rate in at the 
left (at x) is 




The heat transfer rate on the right is 



(2(jc + ik) = Q(x) + 



dQ 
dx 



dx + - 



(2.10) 



Using the conditions on the overall heat flow and the expressions in (2.9) and (2. 10) 



Q{x)-\Q{x) + ^{x)dx+-\ = 0. 
ax 



Taking the limit as dx approaches zero we obtain 
dQix) 



(2.11) 



dx 



= 0, 



(2.12a) 



or 
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(2.12b) 



If A: is constant (i.e. if the properties of the bar are independent of temperature), this reduces to 



ax\ ax 



(2.13a) 



or (using the chain rule) 



dx 



^ ^ A dx dx 



(2.13b) 



Equations (2.13a) or (2.13b) describe the temperature field for quasi-one-dimensional steady state 
(no time dependence) heat transfer. We now apply this to some examples. 

Example 2.1 : Heat transfer through a plane slab 



T= Tj 




Slab 



jc= 



x = L 



->x 



Figure 2.3: Temperature boundary conditions for a slab 



For this configuration, the area is not a function of x, i.e. A = constant. Equation (2. 13) thus became 



d^T 
dx^ 



= 



(2.14) 



Equation (2.14) can be integrated immediately to yield 
dT 



dx 



= a 



(2.15) 
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and 



T = ax + b . 



(2.16) 



Equation (2.16) is an expression for the temperature field where a and b are constants of integration. 
For a second order equation, such as (2.14), we need two boundary conditions to determine a and b. 
One such set of boundary conditions can be the specification of the temperatures at both sides of the 
slab as shown in Figure 2.3, say T (0) = Ti; T (L) = T2. 

The condition T (0) = Ti implies that b = Ti. The condition T2= T (L) implies that T2= aL+ Ti, or 

a = — -. 

L 

With these expressions for a and b the temperature distribution can be written as 




(2.17) 



This linear variation in temperature is shown in Figure 2.4 for a situation in which Ti > T2. 

T 
▲ 




X 



Figure 2.4: Temperature distribution through a slab 



The heat flux q is also of interest. This is given by 

q = -k— = -k^^^^ = constant . (2.18) 
dx L 

How specific do we need to be about when the one-dimensional assumption is valid? Is it 
enough to say that dA/dx is small? (MP HT.2) 

Why is the thermal conductivity of light gases such as helium (monoatomic) or hydrogen 
(diatomic) much higher than heavier gases such as argon (monoatomic) or nitrogen 
(diatomic)? (MP HT.3) 
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2.2 Thermal Resistance Circuits 



There is an electrical analogy with conduction heat transfer that can be exploited in problem 
solving. The analog of Q is current, and the analog of the temperature difference, Ti - T2, is voltage 
difference. From this perspective the slab is a pure resistance to heat transfer and we can define 

T - T 

Q = -^ (2.19) 



where R = L/kA, the thermal resistance. The thermal resistance R increases as L increases, as A 
decreases, and as k decreases. 

The concept of a thermal resistance circuit allows ready analysis of problems such as a composite 
slab (composite planar heat transfer surface), hi the composite slab shown in Figure 2.5, the heat 
flux is constant with x. The resistances are in series and sum to R = Ri + R2. If Iz, is the temperature 
at the left, and Tr is the temperature at the right, the heat transfer rate is given by 




(2.20) 



^v\^^wv 

Ri R2 

Figure 2.5: Heat transfer across a composite slab (series thermal resistance) 



Another example is a wall with a dissimilar material such as a bolt in an insulating layer. In 
this case, the heat transfer resistances are in parallel. Figure 2.6 shows the physical configuration, 
the heat transfer paths and the thermal resistance circuit. 
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Ri 

ANV 



Figure 2.6: Heat transfer for a wall with dissimilar materials (Parallel thermal resistance) 

For this situation, the total heat flux Q is made up of the heat flux in the two parallel paths: 
Q = Qi + Qi with the total resistance given by: 



1 _ 1 ^ 1 

R 7?j /?2 



(2.21) 



More complex configurations can also be examined; for example, a brick wall with insulation 
on both sides. 



Brick 




io°c 



R] R2 R3 

• — ^^/\^-•A^/\^-•-A^/v — • 

Ti T2 Ts T4 



\ Insulation 
0.03 m 



r 



Figure 2.7: Heat transfer through an insulated wall 



The overall thermal resistance is given by 



(2.22) 



Some representative values for the brick and insulation thermal conductivity are: 
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^Wcfc = k2 = 0.7 W/m-K 
hnsulation = kj = ks = 0.07 W/m-K 



Using these values, and noting that ^7 = A2 = A3 = A, we obtain: 



AR, =AR,=^= ^-^^"^ = 0.42 KAV 
^ k, 0.07 W/mK 



AR,=^= ^-^"^ =0.14m^KAV 
^2 0.7 W/mK 



This is a series circuit so 



O , , T.-T. 140K 

— = constant throughout = — = 

A RA 0.98m' KAV 



2 




Figure 2.8: Temperature distribution through an insulated wall 

The temperature is continuous in the wall and the intermediate temperatures can be found 
from applying the resistance equation across each slab, since Q is constant across the slab. For 
example, to fmd T2: 

q='^lZR = 142 W/m' 
R,A 

This yields T; - = 60 K or T2 = 90 °C. 

The same procedure gives T3 = 70 °C. As sketched in Figure 2.8, the larger drop is across the 
insulating layer even though the brick layer is much thicker. 

What do you mean by continuous? (MP HT.4) 

Why is temperature continuous in the composite wall problem? Why is it continuous at the 
interface between two materials? (MP HT.5) 
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why is the temperature gradient dT/dx not continuous? (MP HT.6) 

Why is AT the same for the two elements in a parallel thermal circuit? Doesn't the relative 
area of the bolt to the wood matter? (MP HT.7) 

2. 3 Steady Quasi-One-Dimensional Heat Flow in Non-Planar Geometry 

The quasi one-dimensional equation that has been developed can also be applied to non-planar 
geometries. An important case is a cylindrical shell, a geometry often encountered in situations 
where fluids are pumped and heat is transferred. The configuration is shown in Figure 2.9. 



For a steady axisymmetric configuration, the temperature depends only on a single coordinate (r) 
and Equation (2.12b) can be written as 




control volume 



Figure 2.9: Cylindrical shell geometry notation 




(2.23) 



or, since A = 2jir, 




(2.24) 



The steady-flow energy equation (no flow, no work) tells us that Q.^ = Q^^^ or 



dQ_ 
dr 



= 



(2.25) 



The heat transfer rate per unit length is given by 



Q = -k-2jrr 



dT 
dr 
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Equation (2.24) is a second order differential equation for T. Integrating this equation once gives 



r — = a. (2.26) 
dr 

where a is a constant of integration. Equation (2.26) can be written as 

dT = a— (2.27) 

r 

where both sides of equation (2.27) are exact differentials. It is useful to cast this equation in terms 
of a dimensionless normalized spatial variable so we can deal with quantities of order unity. To do 
this, divide through by the inner radius, rj 

dT = af!^ (2.28) 
[r/rj 

Integrating (2.28) yields 

r = aln^^j+Z7. (2.29) 

To find the constants of integration a and b, boundary conditions are needed. These will be taken to 
be known temperatures T] and T2 at rj and respectively. Applying T= Tj at r = ri gives Ti = b. 
Applying T= T2 at r= r2 yields 

7; =aln^ + 7;, 

or 

T -T 
a = ; ' 

The temperature distribution is thus 

(2.30) 

As said, it is generally useful to put expressions such as (2.30) into non-dimensional and 
normalized form so that we can deal with numbers of order unity (this also helps in checking 
whether results are consistent). If convenient, having an answer that goes to zero at one limit is also 
useful from the perspective of ensuring the answer makes sense. Equation (2.30) can be put in non- 
dimensional form as 
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l^ Mrlr\ (2.31) 
T2-T, Hr^/rJ 

The heat transfer rate, Q , is given by 

• , dT ^ (n-T,)\ ^ k(T,-n) 

Q = -kA — =- Inr^k -M -A— = 2jr 4 

dr ln(r2 1 r^Jr^ ln(r2 1 r^) 

per unit length. The thermal resistance R is given by 

R = (2.32) 

2;rA: 

A 7: -r, 



R 

The cylindrical geometry can be viewed as a limiting case of the planar slab problem. To 

r — r 

make the connection, consider the case when — — - « 1 . From the series expansion for + x) 
we recall that 

2 5 

ln(7+x)« X - y +y+--- (2-33) 

(Look it up, try it numerically, or use the binomial theorem on the series below and integrate term by 

term. = l- jc + jc^+...) 
1 + x 

The logarithms in Equation (2.3 1) can thus be written as 

Inli + ^JlZl) ^ and In^^^^ (2.34) 



in the limit of (r^ - r;) « r;. Using these expressions in equation (2.30) gives 

7^ = (7;-7;)t^+7;. (2.35) 

V 2 ^1 / 

With the substitution of r - r; = x, and - r; = Z we obtain 

r = r, + (r2-rj^ (2.36) 
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which is the same as equation (2.17). The plane slab is thus the limiting case of the cylinder if (r - 
ri) I r« 1, where the heat transfer can be regarded as taking place in (approximately) a planar slab. 

ln(l + x) 

To see when this is appropriate, consider the expansion — -, which is the ratio of heat flux for a 

X 

cylinder and a plane slab. 

Table 2.2: Utility of plane slab approximation 



X 


.1 


.2 


.3 


.4 


.5 


ln(l + jc) 

X 


.95 


.91 


.87 


.84 


.81 



For < 10% error, the ratio of thickness to inner radius should be less than 0.2, and for 20% error, the 
thickness to inner radius should be less than 0.5. 

A second example is the spherical shell with specified temperatures T (rj) = Ti and T (r2) = 
T2, as sketched in Figure 2.10. 




Figure 2.10: Spherical shell 

The area is now^(r) = Am-^ , so the equation for the temperature field is 

dry dr ) 

Integrating equation (2.37) once yields 

dT 2 
— = al r . 
dr 

Integrating again gives 
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(2.37) 



(2.38) 



or, normalizing the spatial variable 



r = 



+ b 



(2.39) 



where a' and b are constants of integration. As before, we specify the temperatures at r = r; and r ■ 
r2. Use of the first boundary condition gives T{r^^= = a +b . Applying the second boundary 
condition gives 

Solving for a ' and b, 

T - T 
a'= ' ' 



\-r.lr^ 

T.T (2.40) 
1 - Tj / 

In non-dimensional form the temperature distribution is thus: 

T.-T l-(r/r) 

= — 7^ — -\ (2.41) 
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3.0 Convective Heat Transfer 



The second type of heat transfer to be examined is convection, where a key problem is 
determining the boundary conditions at a surface exposed to a flowing fluid. An example is the wall 
temperature in a turbine blade because turbine temperatures are critical as far as creep (and thus 
blade) life. A view of the problem is given in Figure 3.1, which shows a cross-sectional view of a 
turbine blade. There are three different types of cooling indicated, all meant to ensure that the metal 
is kept at a temperature much lower than that of the combustor exit flow in which the turbine blade 
operates. In this case, the turbine wall temperature is not known and must be found as part of the 
solution to the problem. 



Figure 3.1: Turbine blade heat transfer configuration 



To find the turbine wall temperature, we need to analyze convective heat transfer, which means we 
need to examine some features of the fluid motion near a surface. The conditions near a surface are 
illustrated schematically in Figure 3.2. 



y 



y 





Velocity 
distribution; 
c = at surface 



//////// / ^ 



c (velocity) 



Figure 3.2: Temperature and velocity distributions near a surface. 
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In a region of thickness 8', there is a thin "film" of slowly moving fluid through which most of the 
temperature difference occurs. Outside this layer, J is roughly uniform (this defines 8'). The heat 
flux can thus be expressed as 



Q__±^ (3.1) 

It cannot be emphasized enough that this is a very crude picture. The general concept, however, is 
correct, in that close to the wall, there is a thin layer in which heat is transferred basically by 
conduction. Outside of this region is high mixing. The difficulty is that the thickness of the layer is 
not a fluid property. It depends on velocity (Reynolds number), structure of the wall surface, 
pressure gradient and Mach number. Generally 6' is not known and needs to be found and it is 
customary to calculate the heat transfer using [kjiuid / 8']. This quantity has the sjonbol h and is 
known as the convective heat transfer coefficient. The units of h are W/m^K. The convective heat 
transfer coefficient is defined by 

^ = ^ = A(r,-r.) (3.2) 

Equation 3.2 is often called Newton 's Law of Cooling . For many situations of practical interest, the 
quantity h is still known mainly through experiments. 

Muc/c/y poinf^ 

How do we know that 5' is not a fluid property? (MP HT.8) 
3. 1 The Reynolds Analogy 

We describe the physical mechanism for the heat transfer coefficient in a turbulent boundary layer 
because most aerospace vehicle applications have turbulent boundary layers. The treatment closely 
follows that in Eckert and Drake (1959). Very near the wall, the fluid motion is smooth and laminar, 
and molecular conduction and shear are important. The shear stress, t, at a plane is given by 

iu — = T (where jU is the djoiamic viscosity), and the heat flux hy q = -k — . The latter is the same 
dy dy 

expression that was used for a solid. The boundary layer is a region in which the velocity is lower 

than the free stream as shown in Figures 3.2 and 3.3. In a turbulent boundary layer, the dominant 

mechanisms of shear stress and heat transfer change in nature as one moves away fi^om the wall. 
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Figure 3.3: Velocity profile near a surface 

As one moves away from the wall (but still in the boundary layer), the flow is turbulent. The fluid 
particles move in random directions and the transfer of momentum and energy is mainly through 
interchange of fluid particles, shown schematically in Figure 3.4. 



m' CpT' 




m'cp T 



Figure 3.4: Momentum and energy exchanges in turbulent flow. 

With reference to Figure 3.4, because of the turbulent velocity field, a fluid mass w' penetrates the 

plane aa per unit time and unit area. In steady flow, the same amount crosses aa from the other side. 
Fluid moving up transports heat m' CpT. Fluid moving down transports m' CpT' downwards. If T > 
T', there is a turbulent downwards heat flow, qj, given by qj = m'cp{T' - T) that results. 

Fluid moving up also has momentum m ' c and fluid moving down has momentum m' c'. The net 
flux of momentum down per unit area and time is therefore m' (c' - c). This net flux of momentum 
per unit area and time is a force per unit area or stress, given by 

tj=m'(c'-c) (3.3) 
Based on these considerations, the relation between heat flux and shear stress at plane aa is 
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4f=T/>(^-£) (3.4) 
or (again approximately) 

^t = '^tCp— (3.5) 

since the locations of planes 1-1 and 2-2 are arbitrary. For the laminar region, the heat flux towards 
the wall is 

k dT 
q = r- — 
IX dc 

k 

The same relationship is applicable in laminar or turbulent flow if — = c or, expressed slightly 

H " 

differently, 

k kIpCp a 

where v is the kinematic viscosity, and a is the thermal diffusivity. 

The quantity (ACplk is known as the Prandtl number (Pr), after the man who first presented the idea of 
the boundary layer and was one of the pioneers of modem fluid mechanics. For gases, Prandtl 
numbers are in fact close to unity and for air Pr = 0. 71 at room temperature. The Prandtl number 
varies little over a wide range of temperatures; approximately 3% from 300-2000 K. 

We want a relation between the values at the wall (at which T= and c= 0) and those in the free 
sfream. To get this, we integrate the expression for dT from the wall to the free stream 

dT= — ^dc (3.6) 

Cp X 

where the relation between heat fransfer and shear sfress has been taken as the same for both the 
laminar and the turbulent portions of the boundary layer. The assumption being made is that the 
mechanisms of heat and momentum transfer are similar. Equation (3.6) can be integrated from the 
wall to the freestream (conditions "at oc")\ 

GC 1 ^ 

S^r- Siiy (3.7) 

W -r^ w 

where - and Cp are assumed constant. 

X 
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Carrying out the integration yields 



T^-T^ = ^— (3.8) 

where is the velocity and Cp is the specific heat. In equation (3.8), is the heat flux to the wall 
and Tw is the shear stress at the wall. The relation between skin fi-iction (shear stress) at the wall and 
heat transfer is thus 

1w _ '^w ^2 9) 



2' 

00 



The quantity — - is known as the skin friction coefficient and is denoted by C/ . The skin 

friction coefficient has been tabulated (or computed) for a large number of situations. If we define a 
non-dimensional quantity 

_ h{T^-T^) _ h 



known as the Stanton Number, we can write an expression for the heat fransfer coefficient, h as 

Cf 

h-'p^CpC^^. (3.10) 

Equation (3.10) provides a useful estimate ofh, or q^, based on knowing the skin friction, or drag. 
The direct relationship between the Stanton Number and the skin friction coefficient is 



St = ^ 



2 



The relation between the heat transfer and the skin friction coefficient 

(Iw" 



is known as the Reynolds analogy between shear stress and heat transfer. The Reynolds analogy is 
extremely usefiil in obtaining a first approximation for heat transfer in situations in which the shear 
stress is "known". 
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An example of the use of the Reynolds analogy is in analysis of a heat exchanger. One type of heat 
exchanger has an array of tubes with one fluid flowing inside and another fluid flowing outside, with 
the objective of transferring heat between them. To begin, we need to examine the flow resistance 
of a tube. For fully developed flow in a tube, it is more appropriate to use an average velocity c and 
a bulk temperature Jg . Thus, an approximate relation for the heat transfer is 



vr- p 



(3.11) 



The fluid resistance (drag) is all due to shear forces and is given by t„. = A where A^, is the tube 
"wetted" area (perimeter x length). The total heat transfer, Q, is q^^.A so that 



Q = Dc, 



T -T 



B 



(3.12) 



The power, P, to drive the flow through a resistance is given by the product of the drag and the 
velocity, Dc , so that 



P 



(3.13) 



The mass flow rate is given hy m = pcA where A is the cross sectional area. For given mass flow 
rate and overall heat transfer rate, the power scales as or as 1/A^, i.e. 



(3.14) 



Equations (3.13) and (3.14) show that to decrease the power dissipated, we need to decrease c, 
which can be accomplished by increasing the cross-sectional area. Two possible heat exchanger 
configurations are sketched in Figure 3.5; the one on the right will have a lower loss. 



heat exchanger 



heat exchanger 



m 





-► diffuser 



high loss 



lower loss 



Figure 3.5: Heat exchanger configurations 
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To recap, there is an approximate relation between skin friction (momentum flux to the wall) and 
heat transfer called the Reynolds analogy that provides a useflil way to estimate heat transfer rates in 
situations in which the skin friction is known. The relation is expressed by 



C 



st = 



2 



(3.15a) 



or 



heat flux to wall momentum flux to wall 



(3.15b) 



convected heat flux convected momentum flux 



or 



"^w 



(3.15c) 



2 



The Reynolds analogy can be used to give information about scaling of various effects as well as 
initial estimates for heat transfer. It is emphasized that it is a usefLil tool based on a hypothesis about 
the mechanism of heat transfer and shear stress and not a physical law. 

What is the "analogy" that we are discussing? Is it that the equations are similar? (MP HT.9) 
In what situations does the Reynolds analogy "not work"? (MP HT.IO) 

3.2 Combined Conduction and Convection 

We can now analyze problems in which both conduction and convection occur, starting with a wall 
cooled by flowing fluid on each side. As discussed, a description of the convective heat transfer can 
be given explicitly as 



This could represent a model of a turbine blade with internal cooling. Figure 3.6 shows the 
configuration. 



Q 

A 



(3.16) 
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Figure 3.6: Wall with convective heat transfer 



The heat transfer in fluid 1 is given by 



Q 

A 



which is the heat transfer per unit area to the fluid. The heat transfer in fluid (2) is similarly given by 



- = hAT2-T^2)- 
A 

Across the wall, we have 

^. = -JJw2-Twl)- 
A L 

The quantity QIA is the same in all of these expressions. Putting them all together to write the 
known overall temperature drop yields a relation between heat transfer and overall temperature drop, 

T2-T1 ■ 

We can define a thermal resistance, R, as before, such that 

R 
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where R is given by 



R = — + — + — . (3.18) 

hjA Ak ^ ^ 

Equation (3.18) is the thermal resistance for a solid wall with convection heat transfer on each side. 

For a turbine blade in a gas turbine engine, cooling is a critical consideration. In terms of Figure 3.6, 
T2 is the combustor exit (turbine inlet) temperature and Ti is the temperature at the compressor exit. 
We wish to find T because this is the highest metal temperature. From (3.17), the wall 
temperature can be written as 

Ah2 R Ah2 

Using the expression for the thermal resistance, the wall temperatures can be expressed in terms of 
heat transfer coefficients and wall properties as 



K2=T2- r. (3.20) 



— ^ + 1 



n 



Equation (3.20) provides some basic design guidelines. The goal is to have a low value of T^^ . 

This means hi/h2 should be large, k should be large (but we may not have much flexibility in choice 
of material) and L should be small. One way to achieve the first of these is to have h2 low (for 
example, to flow cooling air out as in Figure 3. 1 to shield the surface). 

A second example of combined conduction and convection is given by a cylinder exposed to a 
flowing fluid. The geometry is shown in Figure 3.7. 




Figure 3.7: Cylinder in a flowing fluid 
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For the cylinder the heat flux at the outer surface is given by 

q = ^ = liT^-T^ at r = r^ 
A 

The boundary condition at the inner surface could either be a heat flux condition or a temperature 
specification; we use the latter to simplify the algebra. Thus, T = JJ at r = . This is a model for the 
heat transfer in a pipe of radius ri surrounded by insulation of thickness r2-ri. The solution for a 
cylindrical region was given in Section 2.3 as 

r 

T(r) = a In — + b 

n 

Use of the boundary condition T{ri) = Ti yields b = Ti. 

At the interface between the cylinder and the fluid, r = r2, the temperature and the heat flow are 
continuous. (Question: Why is this? How would you argue the point?) 



q = -k— = -k- = h\{aln^-^ + Ti\-Tj. (3.21) 
dr r2 [V r? / J 

\ V 

heat flux just inside cylinder surface heat flux to fluid 

Plugging the form of the temperature distribution in the cylinder into Equation (3.21) yields 

-ci— + hln^A = h{Ti -T^). 
The constant of integration, a, is 



a = 



^2 fj Vj 



and the expression for the temperature is, in normalized non-dimensional form 



ln(r/r;) 



(3.22) 
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The heat flow per unit length, Q, is given by 



Q= Mry-r^)^ (323) 

— +ln(/'y/7) 
hr2 

The units in Equation (3.23) are W / m-s. 

A problem of interest is choosing the thickness of insulation to minimize the heat loss for a fixed 
temperature difference Ti - between the inside of the pipe and the flowing fluid far away from 
the pipe. (Ti - is the driving temperature distribution for the pipe). To understand the behavior of 
the heat transfer we examine the denominator in Equation (3.23) as varies. The thickness of 
insulation that gives maximum heat transfer is given by 

-^[A + ln!l]=0 (3.24) 
dr2\hr2 rj ) 

(Question: How do we know this is a maximum?) 

From Equation (3.24), the value of r2 for maximum g is thus 

maximum heat transfer ' = k/h. (3.25) 

If r2 is less than this, we can add insulation and increase heat loss. To understand why this occurs, 
consider Figure 3.8, which shows a schematic of the thermal resistance and the heat fransfer. As 
increases from a value less than r2 = Mi, two effects take place. First, the thickness of the insulation 
increases, tending to drop the heat transfer because the temperature gradient decreases. Secondly, 
the area of the outside surface of the insulation increases, tending to increase the heat transfer. The 
second of these is (loosely) associated with the k/hr2 term, the first with the \n{r2/ri) term. There are 
thus two competing effects which combine to give a maximum Q ditr2 = k/h. 
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R 




Q 





k 
h 



Figure 3.8: Critical radius of insulation 



In the_expression 



"h.A 



1 



. what is A? (MP HT,1 1 ) 



It seems_that we haye_sjm cpny_ection _a_lpt. Is finding the hea_t_trarisfer cgeffici_enti h^ 

reaUyidifficult? (MP HT]12) 

Whatjdoes the "K" in the contact resi_stance _f grmula stand _f or? 1MP_HT. 13) 

In the_equatipn for the temperature. iA a.cyjinder _(3._22)^ what is "r"_?_CMP_ HT.14) 

3-3 _ Dimensipnless 

Phenomena in fluid flow and heat transfer depend on dimensionless_ parameters. The Mach number 

and the Reynolds number are two you have already seen. These parameters give information as to 
the relevant flow regimes of a given solution. Casting equations in dimensionless form helps show 
the generality of application to a broad class of situations (rather than just one set of dimensional 
parameters). It is generally good practice to use non-dimensional numbers, forms of equations, and 
results presentation whenever possible. The results for heat transfer from the cylinder are already in 
dimensionless form but we can carry the idea even further. For the cylinder: 




(3.26) 



The parameter or — , where Z is a relevant length for the particular problem of interest, is 

k k 

called the Biot number denoted by Bi. In terms of this parameter. 
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The size of the Biot number gives a key to the regimes in which different features are dominant. For 
Bi » 1 the convection heat transfer process offers little resistance to heat transfer. There is thus 
only a small outside (i.e. T (r2) close to r„) compared to the 47 through the solid with a limiting 
behavior of 

T -Tj \nr/rj 
T^-T \m-2/rj 

as Bi goes to infinity. This is much like the situation with an external temperature specified. 

For Bi« 1 the conduction heat transfer process offers little resistance to heat transfer. The 
temperature difference in the body (i.e. fi-om r; to r^) is small compared to the external temperature 
difference, Ti-T^. In this situation, the limiting case is 

^^^ = Bi\n(r/r) 

In this regime there is approximately uniform temperature in the cylinder. The size of the Biot 
number thus indicates the regimes where the different effects become important. 

Figure 3.9 shows the general effect of Biot number on temperature distribution. Figure 3.10 is a plot 
of the temperature distribution in the cylinder for values of 5/ = 0.1, 1.0 and 10.0. 
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Figure 3.9: Effect of the Biot Number [hL/kbody] on the temperature distributions in the solid 
and in the fluid for convective cooling of a body. Note that kbody is the thermal conductivity of 

the body, not of the fluid. 

[adapted from: A Heat Transfer Textbook, John H. Lienhard, Prentice-Hall Publishers, 1980] 
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Configuration Solution 



Figure 3.10: Temperature distribution in a convectively cooled cylinder for different values of 
Biot number, Bi; ¥2/ ri = 2 [from: A Heat Transfer Textbook, John H. Lienhard] 



4.0 Temperature Distributions in the Presence of Heat Sources 

There are a number of situations in which there are sources of heat in the domain of interest. 
Examples are: 

1) Electrical heaters where electrical energy is converted resistively into heat 

2) Nuclear power supplies 

3) Propellants where chemical energy is the source 

These situations can be analyzed by looking at a model problem of a slab with heat sources a (W/m^) 
distributed throughout. We take the outside walls to be at temperature and we will determine the 
maximum internal temperature. 



Slice atx, x+dx 




(a) (b) 

Figure 4.1: Slab with heat sources (a) overall configuration, (b) elementary slice 

With reference to Figure 4.1(b), a steady-state energy balance yields an equation for the heat flux, q . 
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q+ adx - 



dq 
q+ ——dx 
dx 



\ 



= 



(4.1) 



/ 



or 



d q 
dx 



= a. 



(4.2) 



There is a change in heat flux with x due to the presence of the heat sources. The equation for the 
temperature is 



d^T 
dx^ 



+ a/k = 



(4.3) 



Equation (4.3) can be integrated once, 



dT a 

— = X + a 

dx k 



(4.4) 



and again to give 



a 



T = - — ■¥ ax + b 
2k 



(4.5) 



where a and b are constants of integration. The boundary conditions imposed are r(o)= t{l)= . 

oL 

Substituting these into Equation (4.5) gives b = T^ and a = — . The temperature distribution is thus 

2k 



ox a 

T = - + — Lx + 

2k 2k 



(4.6) 



Writing (4.6) in a normalized, non-dimensional fashion gives a form that exhibits in a more useful 
manner the way in which the different parameters enter the problem: 



T -T^r 1 ( X x^\ 



(4.7) 



This distribution is sketched in Figure 4.2. It is sjmimetric about the mid-plane at jc = — , with half 

2 

the energy due to the sources exiting the slab on each side. 
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al^ /k . 



0.125 




X 
L 



Figure 4.2: Temperature distribution for slab with distributed heat sources 



The heat flux at the side of the slab, x= 0, can be found by differentiating the temperature 



distribution and evaluating atx = : -k 



dx 



= .k^(-\=-aL/2. 



x=0 



2K\L 



This is half of the total heat generated within the slab. The magnitude of the heat flux is the same at 
x= L, although the direction is opposite. 

A related problem would be one in which there were heat flux (rather than temperature) boundary 
conditions atx= and x= L,so that T„. is not known. We again determine the maximum 
temperature. Atx = and L, the heat flux and temperature are continuous so 



-k— = h{T-Tj at x = 0,L. (4.8) 

dx ^ ' 

Referring to the temperature distribution of Equation (4.6) gives for the two terms in Equation (4.8), 

(4.9) 
(4.10) 



= ki -— + a\ = (-or + /fa) 



h[T-Tj) = h^^^ + ax + b-T} 



Evaluating (4. 10) at x = and L allows determination of the two constants a and b. This is left as an 
exercise for the reader. 

For an electric heated strip embedded betw^een tw^o layers, w^hat would the temperature 
distribution be if the tw^o side temperatures w^ere not equal? (MP HT.15) 
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5.0 Heat Transfer From a Fin 



Fins are used in a large number of applications to increase the heat transfer from surfaces. 
T5^ically, the fm material has a high thermal conductivity. The fm is exposed to a flowing fluid, 
which cools or heats it, with the high thermal conductivity allowing increased heat being conducted 
from the wall through the fin. The design of cooling fins is encountered in many situations and we 
thus examine heat transfer in a fin as a way of defining some criteria for design. 

A model configuration is shown in Figure 5.1. The fin is of length L. The other parameters of the 
problem are indicated. The fluid has velocity c„ and temperature T^. 




Figure 5.1: Geometry of heat transfer fin 



We assume (using the Reynolds analogy or other approach) that the heat transfer coefficient for the 
fm is known and has the value h. The end of the fin can have a different heat transfer coefficient, 
which we can call hi. 

The approach taken will be quasi-one-dimensional, in that the temperature in the fin will be assumed 
to be a fiinction of x only. This may seem a drastic simplification, and it needs some explanation. 

With a fin cross-section equal to A and a perimeter P, the characteristic dimension in the transverse 
direction i&A/P (For a circular fin, for example, A/P = r / 2). The regime of interest will be taken 

h(A/ P) 

to be that for which the Biot number is much less than unity, Bi = — -« 1 , which is a realistic 

k 

approximation in practice. 

The physical content of this approximation can be seen from the following. Heat transfer per unit 
area out of the fin to the fluid is roughly of magnitude ~ h{Tw- T^) per unit area . The heat transfer 
per unit area within the fln in the transverse direction is (again in the same approximate terms) 
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AlP 



, where Ti is an internal temperature. These two quantities must be of the same 



Al P T - T 

magnitude. If h « 1 , then— ^ — « 1 . In other words, if Bi « 1, there is a much larger 

capability for heat transfer per unit area across the fm than there is between the fm and the fluid, and 
thus little variation in temperature inside the fin in the transverse direction . To emphasize the point, 
consider the limiting case of zero heat transfer to the fluid i.e., an insulated fm. Under these 
conditions, the temperature within the fin would be uniform and equal to the wall temperature. 

If there is little variation in temperature across the fin, an appropriate model is to say that the 
temperature within the fin is a function of x only, T= T(x), and use a quasi-one-dimensional 
approach. To do this, consider an element, dx, of the fin as shown in Figure 5.2. There is heat flow 

of magnitude O at the left-hand side and heat flow out of magnitude O =0 + ^^cfcc at the 

ax 

right hand side. There is also heat transfer around the perimeter on the top, bottom, and sides of the 
fin. From a quasi-one-dimensional point of view, this is a situation similar to that with internal heat 
sources, but here, for a cooling fin, in each elemental slice of thickness dx there is essentially a heat 

sink of magnitude Pdxh [t -T^, where Pdx is the area for heat tiansfer to the fluid. 









Qin ^ 




► Qout 









X Xr^dx 

Figure 5.2:Eleinent of fin showing heat transfer 

The heat balance for the element in Figure 5.2 can be written in terms of the heat flux using Q = qA, 
for a fin of constant area: 



qA = P/i{T- T^)dx + [qA+^ dxJ^ 



(5.1) 



From Equation (5. 1) we obtain 



A^ + P/iT-T^) = 0. 
dx 



(5.2) 



In terms of the temperature distribution, T(x): 



d T Ph 



dx 



2 Ak 



[t-tJ) = o. 



(5.3) 
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The quantity of interest is the temperature difference {T - T^), and we can change variables to put 
Equation (5.3) in terms of this quantity using the substitution 



^(T-T^) = ^. (5.4) 
ax ax 



Equation (5.3) can therefore be written as 



,(T-TJ-^(T-TJ = 0. (5.5) 



dx^ 

Equation (5.5) describes the temperature variation along the fin. It is a second order equation and 
needs two boundary conditions. The first of these is that the temperature at the end of the fin that 
joins the wall is equal to the wall temperature. (Does this sound plausible? Why or why not?) 

(T-TJ,__o = To-T^. (5.6) 

The second boundary condition is at the other end of the fin. We will assume that the heat transfer 
fi-om this end is negligible^ The boundary condition at x = Z is 

The last step is to work in terms of non-dimensional variables to obtain a more compact description. 

T -T ~ 

In this we define — as A J , where the values of A T range fi^om zero to one and ^ = xl L, 

where ^ also ranges over zero to one. The relation between derivatives that is needed to cast the 
equation in terms of ^ is 

d d^ d \ d 



dx dx d^ L d^ 

Equation (5.5) can be written in this dimensionless form as 

^-C^AaT^O. (5.8) 



There is one non-dimensional parameter in Equation (5.8), which we will call m and define by 



' Note: We don't need to make this assumption, and if we were looking at the problem in detail we would solve it 
numerically and not worry about whether an analytic solution existed. In the present case, developing the analytic 
solution is useful in presenting the structure of the solution as well as the numbers, so we resort to the mild fiction 
of no heat transfer at the fm end. We need to assess, after all is said and done, whether this is appropriate or not. 
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^2^2^ hPL^ (5.9) 
kA 

The equation for the temperature distribution we have obtained is 

^^-mZ^Ar = 0. (5.10) 

This second order equation has the solution 

KT = ae'"'^+be-'"'^ . (5.11) 



(Try it and see). The boundary condition at § = is 

AT(0) = a + b = l. (5.12a) 

The boundary condition at ^ = i is that the temperature gradient is zero or 
dAT 



d% 



-(L) = mLae -mLbe =0. (5.12b) 



Solving the two equations given by the boundary conditions for a and b gives an expression for A T 

I X -x\ 

in terms of the hyperbolic cosine or cosh: |^cosh x = j 

coshmZ(l-S) 

AT 5^ '- (5.13) 

cosh mL 

This is the solution to Equation (5.8) for a fm with no heat transfer at the tip. In terms of the actual 
heat transfer parameters it is written as 



coshffy-^Vfc^ 



r-r U z)ijL4 J 

^ ^ —= ^. (5.14) 

cosh . — L 
\\kA ) 
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The amount of heat removed from wall due to the fin, which is the quantity of interest, can be found 
by differentiating the temperature and evaluating the derivative at the wall, x= 0: 




x=0 



(5.15) 



or 



QL 



dDT 



mZsinh mL 



= mL tanh mL 



(5.16) 



kA[To-Tj^ 



1 = 



cosh mL 



Q 



= tanh mL 



(5.16a) 



4kAhP[To-T^) 



The solution is plotted in Figure 5.3, which is taken from the book by Lienhard. Several features of 

the solution should be noted. First, one does not need fins which have a length such that m is much 
greater than 3. Second, the assumption about no heat transfer at the end begins to be inappropriate 
as m gets smaller than 3, so for very short fins the simple expression above would not be a good 
estimate. We will see below how large the error is. 
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Figure 5.3: The temperature distribution, tip temperature, and heat flux in a straight one- 
dimensional fin with the tip insulated. [Adapted from: Lienhard, ^ Heat Transfer Textbook, Prentice- 
Hall publishers] 



d^T d^(T-T^) 

Why did you change the variable and write the derivative — — ^ — in the equation 

dx dx 

for heat transfer in the fin? (MP HT.16) 

What types of devices use heat transfer fins? (MP HT.17) 

Why did the Stegosaurus have cooling fins? Could the stegosaurus have "heating fins"? 
(MP HT.18) 
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6.0 Transient Heat Transfer (Convective Cooling or Heating) 



All the heat transfer problems we have examined have been steady state, but there are often 
circumstances in which the transient response to heat transfer is critical. An example is the heating 
up of gas turbine compressors as they are brought up to speed during take-off. The disks that hold 
the blades are large and take a long time to come to temperature, while the casing is thin and in the 
path of high velocity compressor flow and thus comes to temperature rapidly. The result is that the 
case expands away from the blade tips, sometimes enough to cause serious difficulties with 
aerodynamic performance. 

To introduce the topic as well as to increase familiarity with modeling of heat transfer problems, we 
examine a lumped parameter analysis of an object cooled by a stream. This will allow us to see 
what the relevant non-dimensional parameters are and, at least in a quantitative fashion, how more 
complex heat transfer objects will behave. We want to view the object as a "lump" described by a 
single parameter. We need to determine when this type of analysis would be appropriate. To 
address this, consider the temperature difference Ti - between two locations in the object, as 
shown in Figure 6. 1 . 




object 



Figure 6.1: Temperature variation in an object cooled by a flowing fluid 

If the heat transfer within the body and from the body to the fluid are of the same magnitude, 

h[T,-TJi^^[T,-Tj) (6.1) 

where Z is a relevant length scale, say half the thickness of the object. The ratio of the temperature 
difference is 

Ti-T, hL 

K-T^~ k (6.2) 
If the Biot number is small the ratio of temperature differences described in Equation (6.2) is also 

{t,-t\ , , 




« 1. We can thus say \ Ti-T\«T^- and neglect the temperature non-uniformity 



within the object. 
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The approximation made is to view the object as having a spatially uniform temperature that is a 
flinction of time only. Explicitly, T= T(t). The first law applied to the object is (using the fact that 
for solids Cp = Cv = c) 

Qin'PVc^ (6-3) 



where p is the density of the object and Fis its volume. In terms of heat transferred to the fluid, 

Qout =-pVc^^. The rate of h 
dt 

evolution of the temperature is 



Qout = -pVc — . The rate of heat transfer to the fluid is Ah (T - TJ, so the expression for the time 
dt 



Ah(T-T^) = -pVc^. (6.4) 



dt 



The initial temperature, T(0), is equal to some known value, which we can call T, . Using this, 

, T -T 

Equation (6.4) can be written in terms of a non-dimensional temperature difference 



T-T 

I 00 / 



d(T-TA^hA^lT-T^ (6.5) 



seen before, I ^ + — = 1 which has the solution x = ae "'^ . For the present problem the form is 



At time t= 0, this non-dimensional quantity is equal to one. Equation (6.5) is an equation you have 

dx X 
dt X 

T — -hAt/pVc /:\ 

= ae . (6.6) 

The constant a can be seen to be equal to unity to satisfy the initial condition. This form of equation 

pVc 



implies that the solution has a heat transfer "time constant" given by t = 



hA 



The time constant, t, is in accord with our "intuition"; high density, large volume, or high specific 
heat all tend to increase the time constant, while high heat transfer coefficient and large area will 

tend to decrease the time constant. This is the same form of equation and the same behavior you 
have seen for the R-C circuit, as shown schematically below. The time dependence of the voltage in 

dE E 

the R-C circuit when the switch is opened suddenly is given by the equation — + = . There 

dt RC 

are, in fact, a number of physical processes which have (or can be modeled as having) this type of 
exponentially decaying behavior. 
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C 



R 



Figure 6.2: Voltage change in an R-C circuit 



In equation 0^,^ = p.V.c. 



dT 
dt 



(6.3), what is c? (MP HT.19) 



In the lumped parameter transient heat transfer problem, does a high density "slow down" 
heat transfer? (MP HT.20) 



7.0 Some Considerations in Modeling Complex Physical Processes 

In Sections 5 and 6, a number of assumptions were made about the processes that we were 
attempting to describe. These are all part of the general approach to modeling of physical systems. 
The main idea is that for engineering systems, one almost always cannot compute the process 
exactly, especially for fluid flow problems. At some level of detail, one generally needs to model, 
i.e. to define some plausible behavior for attributes of the system that will not be computed. 
Modeling can span an enormous range from the level of our assumption of uniform temperature 
within the solid object to a complex model for the small scale turbulent eddies in the flow past a 
compressor blade. In carrying out such modeling, it is critical to have a clear idea of just what the 
assumptions really mean, as well as the fidelity that we ascribe to the descriptions of actual physical 
phenomena, and we thus look at the statements we have made in this context. 

hL 

One assertion made was that because — « 1 and on the basis of a heat balance, 

k 



j(Tc-T^)-h(T^-T^ 



QO 



we could assume 



body 



interior 



T ) 



«1 
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Based on this, we said that Jbody is approximately uniform and 7^ « ^ody interior- Another aspect is 

hA 

that setting « 1 erases any geometrical detail of the fin cross section. The only place where P 

and A enter the problem is in a non-dimensional combination. A third assumption, made in the fin 
problem, was that the heat transfer at the far end can be neglected. The solution including this 

effect, where is an axial Biot number is given as Equation (7. 1). If the quantity is small, 

k k 
you can see that Equation (7. 1) reduces to the previous result (5.13). 

r - T; _ cosh otZ(i - + / m)sinh otZ(i - 



and 



cosh mL + (5/3x^1 / mZ) sinh mL 

Q (Bia,iai/mL) + tanh mL 



4kAhP(To-T^) 7 + ^^tanhmZ 

mL 



8.0 Heat Exchangers 

The general function of a heat exchanger is to transfer heat from one fluid to another. The basic 
component of a heat exchanger can be viewed as a tube with one fluid running through it and 
another fluid flowing by on the outside. There are thus three heat transfer operations that need to be 
described: 

1) Convective heat transfer from fluid to the inner wall of the tube 

2) Conductive heat transfer through the tube wall 

3) Convective heat transfer from the outer tube wall to the outside fluid 

Heat exchangers are typically classified according to flow arrangement and type of construction. 
The simplest heat exchanger is one for which the hot and cold fluids move in the same or opposite 
directions in a concentric tube (or double-pipe) construction. In the parallel-flow arrangement of 
Figure 8.1a, the hot and cold fluids enter at the same end, flow in the same direction, and leave at the 
same end. In the counterflow arrangement of Figure 8.1b, the fluids enter at opposite ends, flow in 
opposite directions, and leave at opposite ends. 

Alternatively, the fluids may be in cross flow (perpendicular to each other), as shown by the 
finned and unfinned tubular heat exchangers of Figure 8.2. The two configurations differ according 
to whether the fluid moving over the tubes is unmixed or mixed. In Figure 8.2a, the fluid is said to 
be unmixed because the fins prevent motion in a direction (y) that is transverse to the main- flow 
direction (x). In this case the fluid temperature varies with x and y. In confrast, for the unfinned tube 
bundle of Figure 8.2b, fluid motion, hence mixing, in the transverse direction is possible, and 
temperature variations are primarily in the main-flow direction. Since the tube flow is unmixed. 
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both fluids are unmixed in the finned exchanger, while one fluid is mixed and the other unmixed in 
the unfinned exchanger. 




(a) (b) 
Figure 8.1: Concentric tubes heat exchangers, (a) Parallel flow, (b) Counterflow 




To develop the methodology for heat exchanger analysis and design, we look at the problem of heat 
transfer fi-om a fluid inside a tube to another fluid outside. 
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T2 



Figure 8.3: Geometry for heat transfer between two fluids 

We examined this problem before in Section 3.2 and found that the heat transfer rate per unit length 
is given by 

■ Tjtkij, - r„ ) 

+ m — + 

It is useful to define an overall heat transfer coefficient ho per unit length as 

Q = ljjtrjk,i:r,-T,). (8.2) 

From (8.1) and (8.2) the overall heat transfer coefficient, ho , is 
1 r, , 1 

T ^+fi^-+7r- (^-3) 

We will make use of this in what follows. 

A schematic of a counterflow heat exchanger is shown in Figure 8.4. We wish to know the 
temperature distribution along the tube and the amount of heat transferred. 
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a2 



'bi 



t 

Tal 

Figure 8.4: Counterflow heat exchanger 

To address this we start by considering the general case of axial variation of temperature in a tube 
with wall at uniform temperature To and a fluid flowing inside the tube. 



Ti ►" 















x= 




X 



dx 




heating 

Figure 8.5: Fluid temperature distribution along the tube with uniform wall temperature 



The objective is to find the mean temperature of the fluid at x, T(x), in the case where fluid comes in 
at X = with temperature Ti and leaves atx = L with temperature T2 . The expected distribution for 
heating and cooling are sketched in Figure 8.5. 

For heating {To> T), the heat flow from the pipe wall in a length dx is 

q TcDdx = hnD(TQ - T)dx 
where D is the pipe diameter. The heat given to the fluid (the change in enthalpy) is given by 

PUmCp^^dT = mCpdT 
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where p is the density of the fluid, Um is the mean velocity of the fluid, Cp is the specific heat of the 
fluid and m is the mass flow rate of the fluid. 

Setting the last two expressions equal and integrating from the start of the pipe, we find 

}^=}^^dx. (8 
Carrying out the integration, 

fr^cfi St„-T J T„-T 'Ir,' 



I.e. 



ln|^^]=-^^. (8 



Equation (8.6) can be written as 



T -T 

l^-L^e-^^ (8 



where 



This is temperature distribution along the pipe. The exit temperature at x = Z is 

xhDL 

T - T --■ — 

' \ = e '"'^ (8 



T - T 



The total heat transfer to the wall all along the pipe is 

Q = mc^(Tj-T,). 
From Equation (8.9), 
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h%DL 



mcp = 



In 



T0-T2 



The total rate of heat transfer is therefore 
hitDL[Ti - T2) 



Q = 



ln5^ ' (8.10) 



T2-T0 



or Q = hKDLATj^^ 
where ATlm is the logarithmic mean temperature difference , defined as 



Tj-Ty at; - at; 



A^z^ = -^ = ^T^. (8.11) 
In— ln( 



AT/ 



0-^/ \ATj} 

The concept of a logarithmic mean temperature difference is useful in the analysis of heat 
exchangers. We can define a logarithmic mean temperature difference for a counterflow heat 
exchanger as follows. With reference to Figure 8.4, an overall heat balance between the two 
counterflowing streams is 

Q = ^a(:pa[Tal - Tal) = T^bCp^[T,,2 " ^w)- (8-12) 

From a local heat balance, the heat given up by stream a in length dx IS — Jfla c dT^. (There is a 

negative sign since Ta decreases). The heat taken up by stream bis - mbCpi^dTi,. (There is a 
negative sign because T[, decreases as x increases). The local heat balance is: 

- ffia Cp^dT^ = -f^b '^pb'^^b = gdA= qiiDdx (8.13) 
Solving (8.13) for dTa and dTb, we find: 

dT, = -^; dn = -^ (8.14) 
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d[T^-T,) = d^T = - 



( 



qdA = -\— —\q nDdx 



(8.15) 



where W = mc . Also, q= ho AT where ho is the overall heat transfer coefficient. We can then say: 



= -hnnD\ — — I dx . 



Integrating fi-om x = to x = Z gives 



In 



Ta2-Tbl \^_U^A]__A\ 



'^al '^b2 



(8.16) 



Equation (8.16) can also be written as: 



T -T 
T -T 



(8.17) 



where 



a = hn^DL 



J_ 1_ 



We know that 



Tal - Ta2 = Q/ Wa ^2 - Tm = Q/ Wb 



(8.18) 



Thus 



J_ 1_ 



Solving for the total heat transfer: 



Q = 



1 1 



(8.19) 



(8.20) 



Rearranging (8.16) allows us to express | ^ | in terms of other parameters as 

^bl 
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1^1 

1 1\ \Tal-T,2 



(8.21) 



Substituting (8.21) into (8.20) we obtain a final expression for the total heat transfer for a 
counterflow heat exchanger: 

Q^hoKBL ^^"''^''^'^^"''^'^^ (8.22) 



In 



or 

Q = hoKDLAT^M (8.23) 
8. 1 Efficiency of a Counterflow Heat Exchanger 

Suppose we know the two inlet temperatures Tai, Tbi, and we need to find the outlet temperatures. 
From (8.17) 

Ta2 - Tbl = (Tal - Tb2) e " 

Ta2 - Tal = Tbl - Tal + (Tal - Tb2 ) 6 " 

Rearranging (8.18), 

Tb2 ~ Tbl + -r—\Tal ~ Tal ) 



Thus 



Ta2 - Tal = Tbl - Tal + (Tal - Tbl) - ^[t^j - T^^] e 



or 



fc-^.2)=^fc-^M) (8.24) 
where r| is the efficiency of a counterflow heat exchanger: 
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From (8.24) and (8.25) we can find outlet temperatures Tai and r^.- 



nibCpi, nibCpi, 



We examine three examples. 



i) 



mb > ma Cp, 



AT can approach zero at cold end 



( 



r| -» i as hg, surface area , KDLhg 

r„, -T 



Maximum value of ratio 



= 1 



nia c 

Maximum value of ratio ''^ - 



Tai Tbi mb c 



pb 



ii) 



mb Cpt, < ma Cj 



mbC 



a is negative , r] - — — as [ ] ^ °° (j^^ < JF^ ) 



maC 



pa 



Maximum value of 



Tal Tj,i ^ 



pa 



T -T 

Maximum value of — — = 1 

T -T 



Hi) ma Cp^ = mb Cp^ 

temperature difference remains uniform, r]= 1 
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9.0 Radiation Heat Transfer (Heat transfer by thermal radiation) 



All bodies radiate energy in the form of photons moving in a random direction, with random phase 
and frequency. When radiated photons reach another surface, they may either be absorbed, reflected 
or transmitted. The behavior of a surface with radiation incident upon it can be described by the 
following quantities: 

a = absorptance - fraction of incident radiation absorbed 

p = reflectance - fraction of incident radiation reflected 

T = fransmittance - fraction of incident radiation fransmitted 

Figure 9.1 shows these processes graphically. 



— AAAA^ 

Incident radiation 



— ^VV^ V->-Radiiition absorbed, ^ 



Radiation reflected, P "^-^AAA/ — 



Radiation transmitted, 



Figure 9.1: Radiation Surface Properties 

From energy considerations the three coefficients must sum to unity 
a+ p+ r= \ 



(9.1) 



Reflective energy may be either diffuse or specular (mirror-like). Diffuse reflections are 
independent of the incident radiation angle. For specular reflections, the reflection angle equals the 
angle of incidence. 



9. 1 Ideal Radiators 



An ideal thermal radiator is called a "black body". It has several properties: 

1) It has a = 1, and absorbs all radiation incident on it. 

2) The energy radiated per unit area is = oT^ where o is the Stefan-Boltzmann constant. 



o = 5.67 X 10"' 



W 



m K 



(9.2) 



HT-54 



w 



The units of Ej, are therefore — j. 

m 

The energy of a black body, ii is distributed over a range of wavelengths of radiation. We can 
dEh AEh 

define = « , the energy radiated per unit area for a range of wavelengths of width 4A . 
The behavior of e^, is given in Figure 9.2. 



1 2 3 4 5 6 7 
Figure 9.2: Emissive power of a black body at several temperatures - predicted and observed 
i?iT\x^ = 0.2898 cm K [adapted from: A Heat Transfer Textbook by Lienhard, J.] 

The distribution of e^, varies with temperature. The quantity XT ai the condition where e^^ is a 
maximum is given by (AJ) = 0.2898 cm K. As T increases, the wavelength for maximum energy 

emission shifts to shorter values. The frequency of the radiation,/ is given by /= c/A so high energy 
means short wavelengths and high frequency. 




150 
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Figure 9.3: A cavity with a small hole (approximates a black body) 



A physical realization of a black body is a cavity with a small hole. There are many reflections and 
absorptions. Very few entering photons (light rays) will get out. The inside of the cavity has 
radiation which is homogeneous and isotropic (the same in any direction, uniform everywhere). 

Suppose we put a small black body inside the cavity as seen in Figure 9.4. The cavity and the black 
body are both at the same temperature. 




Figure 9.4: A small black body inside a cavity 

The radiant energy absorbed by the black body per second and per in is ajfJ, where Hh the 
irradiance, the radiant energy falling on any surface inside the cavity. The radiant energy emitted by 
the black body is Eb. Since as = 1 for a black body, H = Eb. The irradiance within a cavity whose 
walls are at temperature Tis therefore equal to the radiant emittance of a black body at the same 
temperature and irradiance is a function of temperature only. 

9.2 Kirchhoffs Law and "Real Bodies " 

Real bodies radiate less effectively than black bodies. The measurement of this is the emittance, e , 
defined by 
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Emittance: 



E ^ radiation from real body at T 
El, <- radiation from black body at T 



Values of emittance vary greatly for different materials. They are near unity for rough surfaces such 
as ceramics or oxidized metals, and roughly 0.02 for polished metals or silvered reflectors. A table 
of emittances for different substances is given at the end of this section as Table 9.1, taken from the 
book by Lienhard. 

The level of the emittance can be related to the absorptance using the following arguments. Suppose 
we have a small non-black body in the cavity. The power absorbed per unit area is equal to oH. 
The power emitted is equal to E. An energy balance gives E = EbS = aH = aEb,. Thus 



= a = £ 



(9.3) 



Equation (9.3), the relation a= e , is known as Kirchhoff s Law. It implies that good radiators are 
good absorbers. It was derived for the case when Tbody = T surroundings (cavity) and is not strictly true 
for all circumstances when the temperature of the body and the cavity are different, but it is true if 
= a ,ex = e, so the absorptance and emittance are not functions of A . This situation describes a 
"gray body". Also, since , are properties of the surface, ax= e^^. 



9. 3 Radiation Heat Transfer Between Planar Surfaces 




Surface 2 



Surface 1 



Figure 9.5: Path of a photon between two gray surfaces 



Consider the two intinite gray surfaces shown in Figure 9.5. We suppose that the surfaces are thick 
enough so that, a + p = \ (no radiation transmitted so transmittance = 0). Consider a photon emitted 
from Surface 1 (remembering that the reflectance p=\ - a): 

Surface 1 emits E^ 
Surface 2 absorbs £"7 « 2 
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Surface 2 reflects (l - ) 

Surface 1 absorbs (l - 

Surface 1 reflects E^ (l - ~ ) 

Surface 2 absorbs (l - cfj - )^2 

Surface 2 reflects E^ (l - X^ ~ X^ ~ ^2 ) 

Surface 1 absorbs E^ (l - «2 Xl - «i X^ ~ ^2 )^i 

The same can be said for a photon emitted from Surface 2: 

Surface 2 emits E2 

Surface 1 absorbs ^'jCCi 

Surface 1 reflects £'2 (l - ) 

Surface 2 absorbs £'2(l-ai)bf2 

Surface 2 reflects £'2 (l - a j Xl - « 2 ) 

etc 



We can add up all the energy Ej absorbed in 1 and all the energy E2 absorbed in 2. In doing 
bookkeeping, it is helpful to define P = (l -Wj Xl ~^2 )• The energy Ej absorbed in 1 is 

E^ (1 - a2 )^i + ^1 (1 - c(2 Jbfi (1 ~ ^2 Xl ~ ^1 )+ • • • 

This is equal to 

E,il-a,}x,(l + 13 + 13' +...). 

However 

— = (l-/3)"' =l + /3 + /3' + .... 
1-/3 

E (\ — cc \x 

We thus observe that the radiation absorbed by surface 1 can be written as — "-^^ _ 

\-P 

E (\—(X \x 

Likewise is the radiation generated at 2 and absorbed there as well. 

Putting this all together we find that 



E2- 



( Ej^lj-aj)a2\ E^a^ 



is absorbed by 1 . The net heat flux from 1 to 2 is 
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^netl to2~ ^1 



1-|3 



1-|3 



1 - (l - Ui - + aia2) 



or 

If = , we would have ^ = 0, so from Equation (9.4) 
If body 2 is black, a2=l, and £'2 = oT*. 



Therefore, again, £1 = for any gray surface (Kirchhoff s Law). 
Using Kirchhoff s Law we find, 

^1 ^2 ~ ^1^2 

or, as the final expression for heat transfer between gray, planar, surfaces: 



0(7-/ -7-/) 



^netl to2 ~ 1 1 • i^-^) 

— + — -1 



82 
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Example 1 : Use of a thermos bottle to reduce heat transfer 



Silvered Walls 



Inside of thermos 
(hot fluid) 



Outside of thermos 
(Cold) 



ei= 62= 0.02 for silvered walls 

Ti = 100 °C = 373 K 
T2 = 20 °C = 293 K 



9 1 to 2 ~ 



/ 4 

o(j; -T, 



1100-420 ^^W 

= 6.9 — - 

99 ^ 



0.02 0.02 



For the same AT, if we had cork insulation with k = 0.04, what thickness would be needed? 

kAT kAT (0.04)(80) 
q = so a thickness L = = — — - = 0.47 m would be needed. The thermos is indeed a 



9 



6.9 



good insulator. 



Example 2 : Temperature measurement error due to radiation heat transfer 



Thermocouples (see Figure 9.6) are commonly used to measure temperature. There can be errors 
due to heat transfer by radiation. Consider a black thermocouple in a chamber with black walls. 
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^ 
^ 
^ 

^ 
^ 

^ 
^ 

^ 
^ 

^ 



metal 1 



A" 

N V 



metal 2 




Measured Voltage 



Note: The measured voltage is related to the difference between Ti and T2 
(the latter is a known temperature). 



Figure 9.6: Thermocouple used to measure temperature 



Suppose the air is at 20 °C, the walls are at 100 °C, and the convective heat transfer coefficient is 



What temperature does the thermocouple read? 

We use a heat (energy) balance on the control surface shown in Figure 9.7. The heat balance states 
that heat convected away is equal to heat radiated into the thermocouple in steady state. 
(Conduction heat transfer along the thermocouple wires is neglected here, although it would be 
included for accurate measurements.) 



Heat in (radiation) 




s 

V 
V 

\ 

s 

s 
\ 

V 

s 

V 

\ 

V 

s 
s 

s 




Figure 9.7: Effect of radiation heat transfer on measured temperature. 
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The heat balance is 



- Kir) = 0^(^wa/ - KJ 



) 



(9.6) 



where A is the area of the thermocouple. 



Substituting the numerical values gives 



15[Tj^- 293) = 5.67x10 



■8 



373' -tJ) 



from which we find Ttc = 5 1 °C = 324 K. The thermocouple thus sees a higher temperature than the 
air. We could reduce this error by shielding the thermocouple as shown in Figure 9.8. 



Which bodies does the radiation heat transfer occur between in the thermocouple?(MP 
HT.21) 

Still muddy about thermocouples. (MP HT.22) 

Why does increasing the local flow velocity decrease the temperature error for the 
thermocouple? (MP HT.23) 

9. 4 Radiation Heat Transfer Between Black Surfaces of Arbitrary Geometry 

In general, for any two objects in space, a given object 1 radiates to object 2, and to other places as 
well, as shown in Figure 9.9. 




Radiation shield 



s 
s 
s 
s 
s 
s 
s 
s 
s 
s 
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s 
s 
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s 



Figure 9.8: Shielding a thermocouple to reduce radiation heat transfer error 
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elsewhere 




Figure 9.9: Radiation between two bodies 



We want a general expression for energy interchange between two surfaces at different temperatures. 
This is given by the radiation shape factor or view factor, Fi .j. For the situation in Figure 9. 10 

F^_2 = fraction of energy leaving 1 which reaches 2 
F2.J = fraction of energy leaving 2 which reaches 1 
F]_2 , F2.J are functions of geometry only 



For body 1 , we know that Eb is the emissive power of a black body, so the energy leaving body 1 is 
Ehi A]. The energy leaving body 1 and arriving (and being absorbed) at body 2 is Ebi A 1 F1.2. The 
energy leaving body 2 and being absorbed at body 1 is Eb2A2 F2-1. The net energy interchange from 
body 1 to body 2 is 



A2.T2 




Body 2 



Body 1 



Figure 9.10: Radiation between two arbitrary surfaces 



Ebi Ai F1.2 - Eb2 A2 F2-1 - Q 



(9.7) 



Suppose both surfaces are at the same temperature so there is no net heat exchange. If so, 



Ebi A] F1.2 -Eb2A2 F2-1 = 0, but also Ebi = Eb2. Thus 



A I Fi^2 - A2 F2-1. 



(9.8) 



HT-63 



Equation (9.7) is the shape factor reciprocity relation. The net heat exchange between the two 
surfaces is 

Example: Concentric cylinders or concentric spheres 




Ti 

Figure 9.11: Radiation heat transfer for concentric cylinders or spheres 

The net heat transfer from surface 1 to surface 2 is 

Qi-2 = AFi-2(Ei,1-Ei,2) ■ 

We know that F1.2 = 1, i.e. that all of the energy emitted by 1 gets to 2. Thus 

Ql-2 = M^bl-Eb2) 

This can be used to find the net heat transfer from 2 to 1 . 

Q2-I = ^2F2-j(Eb2 - Ebi) = AjFi-2[Eb2 - Ebi) =^^(^m " 

View factors for other configurations can be found analytically or numerically. Shape factors are 
given in textbooks and reports (they are tabulated somewhat like Laplace transforms), and examples 
of the analytical forms and numerical values of shape factors for some basic engineering 
configurations are given in Figures 9.12 through 9.15, taken from the book by hicropera and DeWitt. 
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Metals 










Nonmetals 






Surface 


Temperature (" C) 




s 


Surface 


Temperature (° C) 




e 


Aluminum 










Asbestos 


40 


0.93 


- 0.97 


Polished, 98% pure 


200 


- 600 


0.04 


- 0.06 


Brick 








Commercial sheet 


90 




0.09 




Red, Rough 


40 


0.93 




Heavily oxidized 


90 


- 540 


0.2 


- 0.33 


Silica 


980 


0.8 


- 0.85 


Brass 










Fireclay 


9980 


0.75 




Highly polished 


260 




0.03 




Ordinary refractory 


1090 


0.59 




Dull plate 


40 


- 260 


0.22 




Magnesite refractory 


980 


0.38 




Oxidized 


40 


- 260 


0.46 


- 0.56 


White refractory 


1090 


0.29 




Copper 










Carbon 








Highly polished electrolytic 


90 




0.02 




Filament 


1040 - 1430 


0.53 




Slightly polished, to dull 


40 




0.12 


- 0.15 


Lampsoot 


40 


0.95 




Black oxidized 


40 




0.76 




Concrete, rough 


40 


0.94 




Gold: pure, polished 


90 


- 600 


0.02 


- 0.035 


Glass 








Iron and Steel 










Smooth 


40 


0.94 




Mild steel, polished 


150 


- 480 


0.14 


- 0.32 


Quartz glass (2mm) 


260 - 540 


0.96 


- 0.66 


Steel, polished 


40 


- 260 


0.07 


- 0.1 


Pyrex 


260 - 540 


0.94 


- 0.74 


Sheet steel, rolled 


40 




0.66 




Gypsum 


40 


0.8 


- 0.9 


Sheet steel, strong rough oxide 


40 




0.8 




Ice 





0.97 


- 0.98 


Cast iron, oxidized 


40 


- 260 


0.57 


- 0.66 


Limestone 


40 - 260 


0.95 


- 0.83 


Iron, rusted 


40 




0.61 


- 0.85 


Marble 


40 


0.93 


- 0.95 


Wrought iron, smooth 


40 




0.35 




Mica 


40 


0.75 




Wrought iron, dull oxidized 


20 


- 360 


0.94 




Paints 








Stainless, poUshed 


40 




0.07 


- 0.17 


Black gloss 


40 


0.9 




Stainless, after repeated heating 


230 


- 900 


0.5 


- 0.7 


White paint 


40 


0.89 


- 0.97 


Lead 










Lacquer 


40 


0.8 


- 0.95 


PoUshed 


40 


- 260 


0.05 


- 0.08 


Various oil paints 


40 


0.92 


- 0.96 


Oxidized 


40 


- 200 


0.63 




Red lead 


90 


0.93 




Mercury: pure, clean 


40 


- 90 


0.1 


- 0.12 


Paper 








Platinum 










White 


40 


0.95 


- 0.98 


Pure, polished plate 


200 


- 590 


0.05 


- 0.1 


Other colors 


40 


0.92 


- 0.94 


Oxidized at 590''C 


260 


- 590 


0.07 


- 0.11 


Roofing 


40 


0.91 




Drawn wire and strips 


40 


- 1370 


0.04 


- 0.19 


Plaster, rough lime 


40 - 260 


0.92 




Silver 


200 




0.01 


- 0.04 


Quartz 


100 - 1000 


0.89 


- 0.58 


Tin 


40 


- 90 


0.05 




Rubber 


40 


0.86 


- 0.94 


Tungsten 










Snow 


10 - 20 


0.82 




Filament 


540 


- 1090 


0.11 


- 0.16 


Water, thickness > 0.1mm 


40 


0.96 




Filament 


2760 




0.39 




Wood 

Oak, planed 


40 

20 


0.8 
0.9 


- 0.9 



Table 9.1: Total emittances for difTerent surfaces [Adapted from: A Heat Transfer Textbook, 

J. Lienhard] 



Figure 9.12: View Factors for Three - Dimensional Geometries 

Figure 9.13: Fig. 13.4~View factor for aligned parallel rectangles 
Figure 9.14: Fig 13.5~View factor for coaxial parallel disk 

Figure 9.15: Fig 13.6— View factor for perpendicular rectangles with a common edge 
[from: Fundamentals of Heat Transfer, F.P. Incropera and D.P. DeWitt, John Wiley and Sons] 
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Muddiest points on heat transfer 

HT.l How do we quantify the contribution o f each mode o f heat trans fer in a given 
situation? 

Developing the insight necessary to address the important parts of a complex situation 
(such as turbine heat transfer) and downplay (neglect or treat approximately) the other 
aspects is part of having a solid working knowledge of the fundamentals. This is an 
important issue, because otherwise every problem will seem very complex. One way to 
sort out what is important is to make order of magnitude estimates (similar to what we 
did to answer when the one-dimensional heat transfer approximation was appropriate) to 
see whether all three modes have to be considered. Sometimes one can rule out one or 
two modes on the basis of the problem statement. For example if there were a vacuum 
between the two surfaces in the thermos bottle, we would not have to consider 
convection, but often the situation is more subtle. 

HT.2 How specific do we need to be about when the one-dimensional assumption is 
valid? Is it enough to say that dA/dx is small? 

The answer really is "be specific enough to enable one to have faith in the analysis or at 
least some idea of how good the answer is". This is a challenge that comes up a great 
deal. For now, if we say that A is an area defined per unit depth normal to the 
blackboard then saying dA/dx is small, which is a statement involving a non-dimensional 
parameter, is a reasonable criterion. 

HT.3 Why is the thermal conductivity o f light gases such as helium (monoatomic) or 
hydrogen (diatomic) much higher than heavier gases such as argon (monoatomic) 
or nitrogen (diatomic)? 

To answer this, we need some basics of the kinetic theory of gases. A reference for this 
is Castellan, Physical Chemistry, Benjamin/Cummings Publishers, 1986. Two factor 

contribute, the collision cross section and the average molecular velocity. For the gases 
mentioned above the dominant factor appears to be the velocity. The kinetic energy per 
molecule at a given temperature is the same and so the lower the molecular weight the 
higher the average molecular velocity. 

HT.4 What do you mean by continuous? 

The meaning is similar to the definition you have seen in the math subjects. A way to 
state it is that the function at a given location has the same value as we approach the 
location independent of the direction we approach from. To say this in a more physical 
manner, the temperature as a function of x has the same value at x = a, say, whether we 
approach location a fi-om the left or from the right. In terms of what the function looks 
like, it will have no "jumps" or discontinuous (step) changes in value. 



HT.5 Why is temperature continuous in the composite wall problem? Why is it continuous 



at the interface between two materials? 



We can argue this point by supposing 7 were not continuous, i.e., there was a different 
temperature on the two sides of an interface between two materials, with the interface 
truly a sharp plane. If so, there would be a finite temperature difference across an 
infinitesimal distance, leading to a very large (infinite in the limit) heat fiux which would 
immediately erase the temperature difference. This argument could also be applied to 
any location inside a solid of uniform composition, guaranteeing that the temperature is 
continuous in each material. 

HT.6 Why is the temperature gradient dT/dx not continuous? 



As derived in class, across an interface the heat flux is continuous. From the first law, for 
a thin control volume that encloses the interface the net heat flow into the control volume 
is zero. As sketched below, the contribution from the heat fiux at the upper and lower 
ends of the control volume is negligible so the heat fiux in one side must be the same as 
the fiux out of the other. The heat fiux, however, is related to the temperature gradient by 

q = -kVT , or, for one-dimension, q^. = -k(—\ . If the heat fiux is continuous, and if 

" \dxJ 

the thermal conductivity, k, is not the same in the two materials, then dT/dx is not 
continuous. 



Material 
A 



in 



Material 
B 



. Control volume 



^out 



Interface 



HT. 7 Why is AT the same for the two elements in a parallel thermal circuit? Doesn't the 
relatiye area of the bolt to the wood matter? 

In terms of the bolt through the wood wall, the approximation made is that the bolt and 
the wood are both exposed to the same conditions at the two sides of the wall. The 
relative areas of the bolt and the wood indeed do matter. Suppose we consider a square 
meter area of wood without bolts. It has a certain heat resistance. If we now add bolts to 



the wall, the resistance of each bolt is i?^^;^ = 



If there are n bolts, they will be 



Holt^bolt 

in parallel, and the effective area will be n times the area for one bolt. The thermal 



resistance of the wood will be increased by a factor of A^nhout I \,ith ' which is larger 



In summary, the amount by which the heat transfer is increased depends on the fractional 
area with low thermal resistance compared to the fractional area with high thermal 
resistance. 

HT.8 How do we know that 6' is not a fluid property? 



The term — represents the resistance to heat transfer for a unit area. The resistance to 



heat transfer per unit area (1/heat transfer coefficient) can be computed for cases of 
laminar flow, or measured experimentally where we cannot compute it, and it is found 
that for the same state variables, it can have a range of values of several orders of 
magnitude depending on the parameters I described (Reynolds number, surface condition, 
surface shape. . .). Put another way, if the value of the resistance is affected by the surface 
condition (smooth, bumpy, corrugated, etc.) how can the resistance be just a property of 



HT.9 What is the "analogy" that we are discussing? Is it that the equations are similar? 

While the equations are similar, the concepts are deeper than that. The analogy is drawn 
between the heat transfer process (transfer of heat represented by heat flux) and the 
momentum transfer process (transfer of momentum represented by shear stress) 

HT.IO In what situations does the Reynolds analogy "not work"? 

The Reynolds Analogy is just that. It is not a law of nature, but rather a plausible 
hypothesis that allows useful estimation of the heat transfer coefficients in many 
situations in which little or no explicit heat transfer information exists. In the form we 
have derived it, the Rejmolds Analogy is appropriate for use in air, but it is not strictly 

fXC 

applicable if there are pressure gradients, or if the Prandtl number (Pr = — -) is not 

k 

unity. However, the conceptual framework provided by the analogy has been found 
useful enough that the analogy has been extended (in a more complex form, as briefly 
discussed in class) for application to these situations. 




than unity. 




the fluid? 



HT.ll In the expression , what is A? 

h.A 



A is the area normal to the heat flow. For a turbine blade, for example, it would be the 
outer surface area of the blade. 

HT.ll It seems that we have simplified convection a lot. Is finding the heat transfer 
coefficient, h, really difficult? 

We have indeed "simplified convection a lot". We will look at heat transfer by 
convection in more depth in a few lectures, but to answer the question in a few sentences, 
finding the heat transfer coefficient is a difficult problem, because it necessitates 
determining the fluid djmamics; the latter is key to predictions of heat transfer. Even 
with present computational power, calculating the flow around aerospace devices with 
the accuracy needed to be confident about heat transfer coefficients is not by any means a 
"standard" calculation. For some circumstances, it is still beyond the state of the art. We 
will concentrate on describing (i) the basic mechanisms of convective heat transfer and 
(ii) ways of estimating the heat transfer coefficient from known fluid dynamic 
information. 

HT. 13 What does the "K" in the contact resistance formula stand for? 

The definition of the resistance comes out oi Q = — . The units of Q are 

^thermal 

2 2 

Watts/meter , so the units of ^thermal [degrees x meters AVatts]. The K is thus the 
symbol for Kelvin, or degree centigrade. 

HT.14 In the equation for the temperature in a cylinder (3.22), what is "r"? 

The variable r denotes the radial coordinate, in other words the location of the point at 
which we want to know the temperature. 

HT.15 For an electiic heated strip embedded between two layers, what would the 
temperature distribution be if the two side temperatures were not equal? 

If the two temperatures on the outer surfaces of the composite layer were and , the 
heat flux from the two sides of the electrical heater would be different, but the sum of the 
two heat fluxes would still be equal to the heat generated per unit area and unit time. In 
fact heat could be coming into one side of the heater if one of the temperatures were high 
enough. (Sketch this out and prove it to yourself) The temperature distribution would 
be linear from the heater temperature to the two surface temperatures. If this example is 
not clear, please come and see me. 



HT.16 Why did you chanse the variable and write the derivative — ^as — — in 
— dx^ ~ dx^ 

llic equation for heat transfer in the fin? 

The heat balance and derivation of the equation for temperature (5.3) is given in Section 
5.0 of the notes. This is 



d^T Ph 



dx^ Ak 



(r-r^) = o (5.3) 



It is not necessary to change variables, and one could solve (5.3) as is 

d^T Ph 

(or in the form — ^ T = T^). However, since: (1) the reference temperature is T^, 

dx Ak 

and what is of interest is the difference T-T^, (2) making the substitution results in a 
simpler form of the equation to be solved, and (3) the derivative of T-T^ is the same as 
the derivative of T, I found it convenient to put it the equation in the form of Eq. (5.5), 



dx^ Ak 



HT. 1 7 What types of devices use heat transfer fins? 

A number of types of heat exchangers use fins. Examples of the use of tins you may 
have seen are cooling fms on motorcycle engine heads, cooling fms on electric power 
transformers, or cooling fins on air conditioners. 

HT.18 Why did the Stegosaurus have cooling fins? Could the stegosaurus have "heating 
fins"? 

My knowledge of this issue extends only to reading about it in the text by Lienhard (see 
reference in notes). The journal article referenced in that text is: 

Farlow, J. O., Thompson, C. V., and Rosner, D.E., "Plates of the Dinosaur Stegosaurus: 
Forced Convection Heat Loss Fins?, Science, vol. 192, no. 4244, 1976, pp. 1123-1125 
and cover. 

dT 

HT.19 In equation 0,^ = p.V.c. — (6.3), what is c? 
d^ 

c is the specific heat for a unit mass. We don't have to differentiate between the two 
specific heats for a solid because the volume changes are very small, unlike a gas. 

HT.20 In the lumped parameter transient heat transfer problem, does a high density 
"slow down " heat transfer? 



It doesn't. The high density slows down the rate at which the object changes 
temperature; high density means more "heat capacity". 



HT.21 Which bodies does the radiation heat transfer occur between in the thermocouple? 

The radiation heat transfer was between the walls (or more generally the boundaries of 
the duct) and the thermocouple. The boundaries may not be at the same temperature of 
the flowing fluid, for example in a turbine. If the boundaries are not at the same 
temperature as the fluid, and they radiate to the thermocouple, there can be an error in the 
temperature that the thermocouple reads. The discussion was about different possible 
sources of this type of error. 

HT.22 Still muddy about thermocouples . 

I didn't mean to strew confusion about these. As on page 59 of the notes, if we take a 
pair of dissimilar wires (say copper and constantan (an alloy of tin and several other 
metals) or platinum and rhodium) which are joined at both ends and subject the two 
junctions to a temperature difference, it is found that a voltage difference will be created. 
If we know the temperature of one junction (say by use of an ice bath) and we know the 
relation between voltage and temperature difference (these have been measured in detail) 
we can find the temperature of the other junction from measurement of the voltage. The 
assumption is that the temperature of the junction is the temperature that is of interest; the 
possible error in this thinking is the subject of muddy point 1. 

HT.23 Why does increasing the local flow velocity decrease the temperature error for the 
thermocouple? 

The heat transferred by convection is Qconvecuon = thermocouple - ^jiuid )■ If we neglect 

conduction from the wire junction (in other words assume the thermocouple wires are 
thin and long), the heat balance is between convection heat transfer and radiation heat 
transfer. For a fixed O , if the heat transfer increases the temperature difference 

i^convection ^ 

between the thermocouple and the fluid decreases. We have seen, however, that the heat 

transfer may be estimated using the Reynolds analogy (Section 3.1). For fixed skin 
Iriction coefficient the higher the velocity the higher the heat transfer coefficient. 



